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Abstract 

A monomer-dimer reaction lattice model with lateral repulsion among the 
same species is studied using a mean-field analysis and Monte Carlo simulations. 
For weak repulsions, the model exhibits a first-order irreversible phase transition 
between two absorbing states saturated by each different species. Increasing the 
repulsion, a reactive stationary state appears in addition to the saturated states. 
The irreversible phase transitions from the reactive phase to any of the saturated 
states are continuous and belong to the directed percolation universality class. 
However, a different critical behavior is found at the point where the directed per- 
colation phase boundaries meet. The values of the critical exponents calculated 
at the bicritical point are in good agreement with the exponents corresponding 
to the parity-conserving universality class. Since the adsorption-reaction pro- 
cesses does not lead to a non-trivial local parity-conserving dynamics, this result 
confirms that the twofold symmetry between absorbing states plays a relevant 
role in determining the universality class. The value of the exponent 82, which 
characterizes the fluctuations of an interface at the bicritical point, supports 
the Bassler-Brown's conjecture which states that this is a new exponent in the 
parity-conserving universality class. 



1 Introduction 



Nonequilibrium models are relevant to a broad scope of phenomena in physics, chem- 
istry, biophysics, ecology, etc. A relevant feature of the nonequilibrium models which 
exhibit a second-order irreversible phase transition (IPT) to a unique absorbing state 
is that their critical behavior is in the directed percolation (DP) universality class [|l|. 
DP critical behavior is observed over wide-ranging problems emerging from different 
disciplines such as quantum particle physics p|, irreversible catalytic systems 0, 
the contact process 0, branching annihilating random walks with an odd number of 
offspring 0, etc. This fact led Janssen and Grassberger to conjecture that in one- 
component models continuous transitions to a single absorbing state are in the DP 
universality class |p. Motivated by these findings several models with multiple ab- 
sorbing states were proposed but no new universality class was found 0. This proved 
that a greater number of absorbing states is not enough to obtain a critical behavior 
different from DP. Thus the DP universality class is apparently extremely robust. 

In contrast to the well-established DP universality class, only a few exceptions are 
known that do not belong to this class. The known examples are models A and B 



of probabilistic cellular automata ||T0|, [TT|, nonequilibrium kinetic Ising models [|T2 



the interacting monomer-dimer model with infinite repulsion ||T3|, and the branching 
annihilating random walks with an even number of offsprings [l^, A relevant 

feature is shared by all these models: the number of particles (or kinks) is conserved 
modulo 2. That is why this new universality class is sometimes called parity-conserving 
(PC) class. 

Recently, a hierarchy of unidirectionally coupled DP processes has been studied 
||17|| . It has been shown by means of field-theoretic renormalization group techniques 
and Monte Carlo simulations that new values of the exponent (3 corresponding to the 
order parameter arise at the multicritical point while the dynamical critical exponents 
take the same values as the ones corresponding to DP. This new result again poses the 
question of whether a few distinct universality classes are enough to characterize the 
critical behavior of nonequilibrium systems [p!7|]. 

Few models that display a phase transition in the PC class where no explicit conser- 
vation law is present have appeared in the literature. The known examples are the gen- 



eralized Domany-Kinzel cellular automata |18], the three-species monomer-monomer 



model [T^l and the monomer- monomer surface reaction model 71] . Since these 
models do not explicitly conserve any quantity modulo 2 they show that rather than 
parity conservation the symmetry among absorbing states is the origin for the emer- 
gence of a different class [|I^, . 

In this work we study the behavior of the one- dimensional monomer-dimer surface 
reaction model with lateral repulsion by means of a mean-field (MF) analysis and Monte 
Carlo simulations. This model is an extension in one dimension of the well-known Ziff, 
Gulari, and Barshad surface-reaction model [H]. The model was first proposed by Kim 



and Park [|T^ and studied in the case of infinitely strong repulsions. For finite repulsions 
the model has a rich critical behavior displaying first- and second-order IPT. 

The manuscript is organized as follows. In section II we begin with a brief descrip- 
tion of the model and show the phase diagram obtained by simulations. In section III 
we present the MF analysis and compare MF results with simulation results. Section 
IV contains a detailed analysis of the critical behavior of the model performed by means 
of Monte Carlo simulations. In the last section we summarize our results. 

2 The Model 



The model we study in this paper was first introduced by Kim and Park [|T3]. The 
monomer-dimer reaction model with lateral interactions can be defined as follows: A 
monomer A or a dimer B2 adsorb at the vacant sites of a one-dimensional lattice with 
probabilities p and q respectively, where p + q = 1. Each adsorption attempt begins 
by selecting one site of the lattice at random and if that site is occupied the trial 
ends. Otherwise, if the site is empty, an A {B2) species is selected with probability p 
(g), respectively. In order to introduce lateral interactions between the same kind of 
particles, it is assumed that the adsorption probability of the selected species depends 
on the configuration of the adsorbed particles on the nearest neighboring (NN) sites of 
the selected one. Then the adsorption probability Pa can be written as 



P 



A 



'p if ^ NN A 

p(l - ta) if 3 NN A, 



where < < 1. For the adsorption of B2 one has first to select at random a NN site 
of the empty one and if that site is occupied the trial ends because the dimer can not 
be deposited on the lattice. Otherwise, if the site is empty, the adsorption probability 



is given by 

= 'f^'^''^ (2) 

' \q{l-rB) if 3 NNB. ^ ^ 

where < < 1. In this way, the parameters va and tb can be interpreted as the 
repulsive interactions among similar species. Unlike species on adjacent sites react 
immediately and leave the lattice, leading to a process limited only by adsorption. 



In this work we study the case = = r. Previous studies of this model [13, 22 



focused on various aspects of the critical behavior for infinite repulsive lateral interac- 
tions, i. e. = = 1. In this case, there are two equivalent absorbing states whose 
configurations are given only by monomers occupying the odd or even-numbered lat- 
tice sites respectively. The critical behavior of this model was found to be in the PC 
universality class. 

Figure |I| shows the phase diagram of the system obtained by both static and dynamical 
Monte Carlo simulations. There are two different absorbing states characterized by the 
lattice saturated by A or 5 species respectively. For weak repulsions the system dis- 
plays first order IPT between the saturated states, and no reactive phase is observed. 
The first order critical points have been calculated by means of static Monte Carlo sim- 
ulations since correlations at first order IPT are short ranged. Increasing the repulsion 
a reactive window appears whose edges are second order critical lines that separate 
this state from the absorbing states. The second order IPT have been calculated using 
dynamic Monte Carlo simulations since fluctuations in this case are imposing. 



For r = 1 (infinite repulsion) one can identify since the very beginning the pres- 
ence of two equivalent absorbing states. Although one can not find a priori a situation 
like this for < r < 1, a phase point where both absorbing states are statistically 
equivalent can exist. This happens exactly where both second order lines meet, i. 
e. at the bicritical point. A similar phase diagram was also found in the interacting 
monomer- monomer reaction model . Very recently, a careful study at the bicritical 
point led to the conclusion that this point is in the PC universality class Given 
the symmetry of the adsorption process in the interacting monomer-monomer reaction 
model the bicritical point was found on the line p = q = 0.5. Since the model we study 
in this work does not display this symmetry, the localization of the bicritical point is 
more difficult. 



3 Mean-Field Theory 



In order to obtain a qualitative understanding of the behavior of the model a mean-field 



analysis is performed. We consider mean-field approximations that study the time 
evolution of blocks up to three sites, neglecting higher-order correlations. 

We start considering a one-dimensional system of size L at time t. Each site can be 
only at three different states, namely. A, B, oiV corresponding to a site occupied by a 
particle of type A, B or empty, respectively. In order to calculate MF approximations 
taking into account correlations up to three lattice sites, we write down the following 
MF rate equations 



^ = PPvvv +p{l-r) {pAVA + Paw + Pvva) - g(l - r) (pbvva + Pavvb) 

—qC^PAvvA + Pavvv + Pvvva), (3) 



^ = qC^Pvvvv + Pavvv + Pvvva) + g(l - r) {pbvvv + Pvvvb + Pbvv + Pvvs) 
- r){pAVB + Pbva) - p{pBVB + pBvv + Pvvb), (4) 



= p(l - r){pvvA + Paw + ^pava) - g(l - r){pAAVVB + Pbvvaa) 

dt 

—q{Paavva + Pavvaa + Paavvv + Pvvvaa), (5) 



dpBB 

dt 



qpvvvv - p(l - r){pBBVA + Pavbb) 

+q{l - r){pBvvB + Pbwv + Pvvvb + Pbvv + Pvvb) 

11 

—PyPVVBB + PBBVV + -PBBVB + -PBVBB), (6) 



dpAAA 



dt 



= p(l - r) {pvvAA + Paavv + Pava + Paava + Pavaa) 
-q{l - r){pBvvAAA + Paaavvb) 
~(1{Paaavva + Paw AAA + Pvvvaaa + Paaavvv)^ 



(7) 



BBB 



dt 



dp_ 



'AVA 



dt 



q{l — r){2pBvvB + Pvvbb + Pbbvv + Pbvvv + Pvvvb) 
- r){pBBBVA + Pavbbb) 

—p{pvVBBB + PBBBVV + -PSVBBB + -PBBBVs), (8) 



PiPAVVV + PVVVa) + 2p(l - r)pAVVA - p(l - r)pAVA 

-q{l - r){pBvvAVA + Pavavvb) 

—QiPAVAvvA + Paw AVA + Pvvvava + Pavavvv), (9) 



dp 



BVB 



dt 



Q{PbVVVV + PVVVVB + PAVVVB + PBFVyA) 

+2g(l - r) Pbvvv B - p(l - r){pBVBVA + PAvsyi?) 
"PlPyvBVB + Pbvbvv + Pbvbvb + Pbvb), (10) 



and 



^^^^-^ = PPvvvB + qipvvvvA + Pavvva) + p(l - r)pAyyi? 

+g(l - r)pAvvvB - p{pAVBVv + 2P^v'ByB) - p(l - ?^)pyiyByyi 
-^(PVFyAVB + PAVVAVb) - g(l - r)pBVVAVB - p(l - ?')PAyi?, (11) 

where Pi^i^is... is the density of the block 11121^... and we have used the equations that 
hnk the density of an n-block and the density of an + 1-block 

Piii2...i„ ~ ^ ] Phi2---inin+i ~ ^ ] Pioiii2...i„ (12) 
in+i io 

where AB pairs are not allowed since they immediately reacts. The processes consid- 
ered to obtain equations (^) to (p!T| ) are listed in table 1. Furthermore, we obviously 
have the following constraint 



PA + Pb + PV = 1- 



(13) 



Table 1: Probabilities for the allowed kinetic processes. 
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In the simple MF analysis we just neglect correlations among sites, i. e. Pi^i2i3...im ~ 
PiiPi2Pi3---Pim- Within this approximation, equations @, @, and (|13D, comprise a close 
set of equations. Then a solution for pA, Pb, and pv can be obtained. However, more 
equations are needed in order to obtain results for higher order MF analysis. In the 
pair MF approach we approximate the density of blocks longer than two sites in the 
following way 

m— 1 

pi 



Pil 



Piii2 \_ \_ 
i=2 



Pu 



(14) 



Using equations (0)-®, ([T3|) and writing down the equations (|T2D which relate single- 
site densities with pair densities a solution in the pair-MF approximation can be ob- 
tained. Proceeding in a similar way, within the three-sites MF approximation density 
of blocks longer than three sites are replaced by 



Piii2i 



m-2 
i=2 



Pijij+iij+2 
Pijij+i 



(15) 



Due to the immediate reaction of AB pairs, the stationary densities of triples AVB and 
BVA are equal. Then, by means of equations (|^) to and (0), and considering 
all the relations between single-site densities, density of pairs and density of triples, 
the three-sites MF approximation can be solved. It should be pointed out that due 
to the complexity of the set of equations all the MF approximations have been solved 
numerically. 

Figure |^ shows plots of the densities of A and B in the stationary state versus p and q 
respectively for two different values of the repulsion r obtained by static Monte Carlo 
simulations, simple, pair, and three-sites MF analysis. For weak repulsion (r = 0.5) 



the sharp jump observed in both the densities of A and B is the signature of first order 
IPT. Since correlations at a first order IPT are finite we expect the MF approaches to 
give good results. Although the simple MF approximation is quite poor, results quickly 
improve for the pair and three-sites approaches. For higher values of the repulsion r the 
sharp variation in the densities is no longer present but a smooth transition is observed. 
In general, MF theories do not give good results near second order continuous IPT since 
second order IPT are governed by fluctuations. However, as it is observed in figures 2c 
and 2d, we still obtain a fairly good agreement between the three-sites MF approach 
and simulation results. 

Figure § shows plots of the phase boundary for both the A and B saturated phases 
obtained by simulations, pair, and three-sites MF approximations. Within the simple 
MF analysis the phase boundary for the A {B) saturated phase always occurs at p = 1 
(g = 1), that is why it is not included in the figure. For weak repulsions, the pair MF 
approach gives better results for both the A and B phase boundaries than the three-sites 
MF analysis. However, for stronger repulsions, correlations become more important 
and the three-sites MF approximation leads to better results. It should be pointed out 
that no bicritical point can be obtained from the MF approximations considered in 
this work. However, it is observed in figure ^ that phase boundary curves qualitatively 
resemble the actual phase diagram. The closest points in the phase diagram obtained 
within the three-site MF approach are (rc ~ 0.6, ~ 0.348), (rc ~ 0.6, ~ 0.305) 
which are good approximations for the actual bicritical point [vc = 0.559, pc = 0.35) 
(see next section). In the MF treatment of the interacting monomer- monomer reaction 
model ||2T| a bicritical point can always be found given the symmetry of the adsorption 
process which is reflected in the MF rate equations. However, no good approximation 
for this point was obtained up to the three-sites MF analysis. 

The monomer-dimer reaction model with lateral interactions displays a feature that 
is not present in the interacting monomer- monomer reaction model. In fact, for q = 1, 
the stationary density is always less than one in spite of the value of the repulsion r. 
However, is a function of r since local configurations like BBVBB are more likely 
to occur when the repulsion is increased. Then, for g = 1, we have a one-dimensional 
random dimer filling problem with lateral interactions. In this case p_B is commonly 
called jamming coverage which we denote as Qj{r). 

Figure ^ shows a plot of Oj(r) versus r obtained by simulations, pair and three-sites 



MF approximations. It should be noted that a MF analysis of the jamming coverage 
requires at least to take into account correlations up to pairs of sites. For r = both 
MF analysis predict 6j(0) = 0.8646..., reproducing the value derived long ago by Flory 
P^ . As it has also been observed for the phase boundary curves, the three-sites (pair) 
MF approximation gives better results for strong (weak) repulsions. 

We did not calculate higher order MF approximations since the algebra becomes 
much more complicated and the MF approaches presented here provide a fairly good 
qualitative understanding of the model. 

4 Simulations Results 

As mentioned in the last section, static Monte Carlo simulations are suitable to obtain 
the coordinates of the first-order transition points. However, second-order IPT are 
dominated by fluctuations, so in a finite system and close to the critical point, the 
stationary states of the reactive phase can irreversibly evolve into the saturated state 
(absorbing state). Due to this circumstance, the precise determination of both critical 
points and critical exponents is rather difficult. However, this shortcoming can be 
avoided by performing an epidemic analysis |T[]. Within this context, the epidemic 
analysis is usually called "defect dynamics" simulations. For this purpose one starts, 
at t = 0, with a small block of vacant sites in an otherwise saturated lattice, i.e. a 
configuration close to one of the absorbing states. Then, the time evolution of this 
block is analyzed by measuring the following properties: (i) The average number of 
vacant sites at time t, N{t), (ii) the survival probability of the block at time t, P{t), 
and (iii) the average distance over which the block has spread at time t, R{t). Finite- 
size effects are absent because the system is taken large enough to avoid the presence 
of vacant sites at the boundaries. For this purpose a lattice of 10'^ cells is enough. 
Averages are taken over 10^ different samples. Near the critical point, the number of 
vacant sites is often very small. Then, we improve the efficiency of the algorithm by 
keeping a list of vacant sites. Time is incremented by l/n{t), where n{t) is the number 
of vacant sites at time t. Time evolution of blocks are monitored up to t = 10^. At 
criticality, the following scaling behavior holds: 

N{t) oc t\ (16) 
P{t) oc t-\ (17) 



and 



R{t) oc t'/^ 



(18) 



where 6, rj and z are dynamic exponents. 

At the bicritical point, it is useful to perform another epidemic analysis often called 
"interface dynamics" simulations. In this case one starts at t = with a minimum 
width interface between two saturated phases. Since AB pairs immediately react, there 
must always be at least a vacancy between two saturated phases. Then, the epidemic 
always survives and consequently 6 = 0. In addition, a second type of "interface 
dynamics" simulation can also be performed. In this case the simulation finishes when 
the interface has collapsed back to its initial width [|1^, At the bicritical point 
equations (|16|), (p!7|), and (|1^) hold but P(t) must be interpreted as the probability that 
the interface has not returned to its minimum width. Interface dynamics simulations 
give us information about the competitive growth of different domains. 
Figure |^ shows log-log plots of N(t), P(t), and R{t) for different values of p close to 
the phase boundary between the i?-saturated and the reactive phase keeping r = 0.9 
constant. The straight line obtained for the three quantities mentioned above at g = 
0.835 is the signature of critical behavior, while slight upward (downward) deviations 
for q = 0.8325 (g = 0.8375) indicate supercritical (sub-critical) behavior, respectively. 
In this way we have determined the critical points along both second order phase 
boundary curves. The analysis at the bicritical point is discussed later. The spreading 
or epidemic analysis is a powerful method since the error bars for the calculations of 
the critical points are on the third digit. The critical exponents obtained at various 
critical points along the phase boundary curves have the same values (within error 
bars) and are in good agreement with the dynamic critical exponents corresponding to 
DP which are the following: 



By drawing the second order phase boundary curves one gets the first clue to the 
position of the bicritical point. The main property of the system at the bicritical point 
is the symmetry of both saturated phases. This means that at the bicritical point both 
absorbing states are statistically equivalent. Then, defect dynamics simulations at the 
bicritical point should give identical results no matter what saturated phase is used to 



6 ^ 0.162, 7] = 0.317, z ^ 1.282 



(19) 



start the simulation. 

Figure ^ shows log-log plots of N{t), P{t), and R{t) obtained from defect dynamics 
simulations started using different absorbing states at the bicritical point. It is observed 
that the critical behavior up to time t ^ 10^ is governed by the same exponents. Then, 
our best estimation of the value of the bicritical point is {r^ = 0.559, = 0.35). It 
should be pointed out that in the study of the interacting monomer-monomer reaction 
model ||2T| the localization of the bicritical point was easier because of the symmetry 



of the adsorption-reaction processes. 

We obtain the following values for the dynamical critical exponents at the bicritical 
point 

S = 0.2910 ± 0.0002, ri = 0.0034 ± 0.0003, z = 1.147 ± 0.0004 (20) 

It should be remarked that the error bars merely indicate the statistical errors obtained 
from regressions. 

The exponent (3 which characterizes the critical behavior of the order parameter 
has been obtained directly in static simulations. For the present model a good choice 
for the order parameter is the average density of empty sites pv- The behavior of the 
order parameter close to a critical point is given by 

Pv ~ |A|^ (21) 

where A is the distance from a point within the reactive phase to the bicritical point. 
Figure |] shows a plot of pv versus A where the points within the reactive phase 
belong to a straight line that bisects the reactive window. Although static simula- 
tions are known to be quite inaccurate to determine (3, we found the reasonable value 
P = 0.88 ±0.005. 

The set of critical exponents calculated at the bicritical point are in good agreement 
with the ones corresponding to the PC universality class. For defect dynamics sim- 
ulations at r > Tc and close to the bicritical point a crossover is observed where the 
critical behavior at short times (long times) is governed by the PC (DP) universality 
class. 

In the following we present the results of the " interface dynamics" simulations. Figure 
^ shows plots of the number of vacancies in the interface N{t) and the average size 
of the interface R{t) obtained using the first type of "interface dynamics" simulation 



(where S = 0). We average over 10^ independent samples and we find the following 
value for the dynamical critical exponents: 



r/i = 0.2840 ± 0.0002, = 1.1506 ± 0.0004 (22) 

This values are in good agreement with the corresponding critical exponents obtained 
in other models in the PC universality class |]13|, 0, |T5|, |18|, |l^. It should be noted 
that the values of rji and zi are similar to the ones obtained from "defect dynamics" 
simulations if we only average over surviving runs. Then the first type of "interface 
dynamics" simulation yields no new information about the dynamics of the interface 
and reflects the full equivalence of the absorbing states at the bicritical point. 

Figure |I0| shows plots of N{t), P{t), and R(t) obtained using the second kind of 
"interface dynamics" simulation. We average over 9 x 10^ independent samples and 
we find the following values for the critical exponents: 

52 = 0.71 63 ± 0.0003, r]2 = -0.4277 ± 0.0004, Z2 = 1.160 ± 0.0006 (23) 

We first observe that the values of the exponents zi and Z2 are similar. The exponent 
?7i must be compared with 62 + ri2 which governs the time evolution of the vacant sites 
averaged only over the surviving runs. Both exponents are equal within error bars. 
This is the signature of the universal behavior of the critical spreading of the interface. 
The value of the new independent dynamical exponent 62 supports the Bassler-Browne 



conjecture [1^ which states that 62 is a new universal exponent within the PC class. 



Recently Park and Park |2^ have found that the interacting monomer-dimer reaction 
model with infinite repulsive interaction (r = 1) supports a kink representation where 
the total number of kinks is conserved modulo 2. By including a parity-conserving 
symmetry breaking field that favors one of the absorbing states, the authors showed 
that the critical behavior of the model changes from PC to DP. Then, they concluded 
that the conservation of the number of kinks modulo 2 is not the reason for observing 
a universality class different from DP but the symmetry of the absorbing states. 

In contrast to the case with infinite repulsion, the present variant of the model 
involve adsorption-reaction processes in which parity is not conserved explicitly. Al- 
though by their own nature the number of domain walls are conserved modulo 2, it 
should be remarked that a non-trivial PC dynamics requires the creation of at least 
three kinks per step which is not possible in our model. This corroborates that the 



symmetry among absorbing states is indeed the only essential property of models in 
the PC class 0, 0, 0. 

It is interesting to discuss the relation between the present model and the branching 
and annihilating walks models. Static Monte Carlo simulations close to the bicritical 
point reveal that stationary configurations are formed by large clusters of different 
species that survive a long time. Then, it is clear that the dynamics of domain walls 
between clusters of different species will be relevant at the bicritical point. Following 



the ideas presented in reference |]T8[, at the bicritical point "walkers" can be identified 
with domain walls between clusters of different species. It should be pointed out that 
we have to generalize the concept of walkers. In fact, the "walkers" defined above are 
extended objects of fiuctuating width. It should be noticed that at the bicritical point 
the typical width of a walker is only few lattice sites. Considering the time evolution 
of this "walkers" in the long-time regime and for large lattices, an effective parity- 
conserving dynamics may be restored. Concerning the dynamics of these "walkers", it 
is possible to identify an inactive and an active phase. In fact, the line defined by the 
first-order phase transition points, corresponds to the inactive phase for the dynamics 
of the "walkers", which ends at the bicritical point. A continuation of this line from the 
bicritical point through the reactive phase, would correspond to the active phase for the 
dynamics (see figure 1). In other words, the dynamics of this walkers is defined within 
a subspace of the phase diagram in which the statistical weight of the kinetic processes 
involving the species A and B is the same. Recently, Cardy and Tauber have developed 
a systematic theory for the branching and annihilating random walkers p5|]. The 
authors have shown that in one dimension, fiuctuation effects lead to the emergence of 
a non-trivial inactive phase for values of the branching rate < a < ac, and a dynamic 
phase transition at o"c > 0, in contrast to the mean-field result Uc = 0. However, 
in two dimensions the theory predicts that fiuctuation effects generate logarithmic 
corrections and the critical branching rate takes the value predicted by the bare mean- 
field theory (ctc = 0). Then, the identification of an inactive phase for < r < Vc and 
a dynamic phase transition at Tc > in the present model, is in complete agreement 
with the theoretical results. Recent results for a two-dimensional monomer-monomer 



surface-reaction model with repulsive lateral interactions pBj] indicate that the critical 
branching rate takes the value Tc = (mean-field value), also supporting the theoretical 
results. 



The same ideas can also be applied to other models in the PC universality class 
that violate local parity conservation []Ty|, pT| . 



5 Conclusions 

We have studied a monomer-dimer surface reaction lattice model with lateral repulsion 
among same species using a mean-field analysis and Monte Carlo simulations. For weak 
repulsions the model exhibits first-order IPT between two phases saturated by different 
species. Increasing the interaction a reactive window appears whose edges are second- 
order critical lines that separate this state from the absorbing states. 

We have considered MF approximations that take into account correlations up to 
three lattice sites. The stationary density of the species, the phase boundary curves, 
and the jamming coverage have been studied within the MF approaches and the results 
compared with simulations. For weak (strong) repulsion the pair (three-sites) MF 
approximation has led us to better results than the three-sites (pair) MF approximation 
for all the quantities mentioned above. No bicritical point has been obtained from 
the MF approximations considered in this work. However, phase boundary curves 
calculated within MF approximations qualitatively resemble the actual phase diagram. 
Within the three-sites MF approach phase boundary curves come up closest at (r^ ~ 
0.6, ^ 0.348), (r, ^ 0.6, ^ 0.305) which are good estimations for the actual 
bicritical point (r^ = 0.559, Pc = 0.35). 

The critical behavior of the model has been studied using both static and dynamical 
Monte Carlo simulations. IPT between the reactive phase and any of the saturated 
phases are in the DP universality class. However, it has been found a critical behavior 
different from DP at the point where both DP curves meet. The critical exponents 
calculated at the bicritical point are in good agreement with the ones corresponding 
to the PC class. Since this model does not conserve explicitly any quantity modulo 
2, it corroborates that the twofold symmetry in the absorbing states is the reason 
for obtaining a critical behavior different from DP. It is relevant to mention that the 
authors of the theory for the branching and annihilating random walkers P31 have 



identified a formal permutation symmetry at the Hamiltonian level in the case of an 
even number of offsprings. This fact again indicates that symmetry among absorbing 
states plays a relevant role in determining the universality class. 

The value of the critical exponent 82 = 0.71 corresponding to the probability P{t) 



that the interface has not collapsed back to its minimum width is in good agreement 
with the values obtained in other models [11151 1211, 17^1 • This corroborates the Bassler- 



Brown conjecture which states that the interfacial fluctuations are an additional 
universal characteristic of models in the PC class. 
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Figure 1: Phase diagram of the model obtained by means of Monte Carlo simulations. 

Figure 2: Plots of and ps versus p and q respectively, obtained using static Monte 
Carlo simulations, simple, pair, and three-sites MF approximations, a) p^ versus p for 
r = 0.5, b) pb versus g for r = 0.5, c) pA versus p for r = 0.9, and d) pb versus q for 
r = 0.9. 

Figure 3: Phase boundary curves obtained by means of Monte Carlo simulations, 
pair, and three-sites MF approximations, a) Boundary of the 74-saturated phase , b) 
Boundary of the S-saturated phase 

Figure 4: Phase boundary curves obtained within the pair, and three-sites MF approx- 
imations, a) Pair MF approximation b) Three-sites MF approximation. 

Figure 5: Plot of the jamming coverage ©j(r) versus r obtained using Monte Carlo 
simulations, pair, and three-sites MF approximations. 

Figure 6: Log-log plots obtained using defect dynamics simulations keeping r — 0.9 
constant, a) N{t) versus t b) P{t) versus t c) R{t) versus t. 

Figure 7: Log-log plots obtained using defect dynamics simulations initiated using 
different absorbing states at the bicritical point, a) N{t) versus t b) P{t) versus t c) 
R{t) versus t. 

Figure 8: Plot of the density of vacant sites pv versus the distance |A| to the bicritical 
point. A straight line of slope (3 — 0.90 has been included. 

Figure 9: Log-log plots obtained using the first type of interface dynamics simulations 
{6 = 0) for the following values of the parameters: (r = 0.57, p = 0.3484) (upper 
curve), (r = 0.559, p = 0.35) (middle curve) (r = 0.559, p = 0.355) (lower curve), a) 
N{t) versus t b) R(t) versus t. 

Figure 10: Log-log plots obtained using the second type of interface dynamics simula- 
tions for the following values of the parameters: (r = 0.57, p = 0.3484) (upper curve), 
(r = 0.559, p = 0.35) (middle curve) (r = 0.559, p = 0.355) (lower curve), a) N{t) 
versus t b) P{t) versus t c) R{t) versus t. The inset amphfies the last decade in order 
to distinguish between the three curves since R{t) is less sensitive to changes in the 
parameters. 
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